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Abstract 

We study matrix models as a new approach to formulate massless higher spin gauge 
field theory. As a first step in this direction, we show that the free equation of motion 
of bosonic massless higher spin gauge fields can be derived from that of a matrix model. 
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1 Introduction 

It has been known that there are problems in the construction of consistent interactions for 
massless higher spin gauge fields, though there are physically acceptable free field Lagrangian 
for them. Free field Lagrangian for higher spin fields, in which these particles are expressed by 
totally symmetric tensor or tensor-spinor fields, was originally derived by Fronsdal for bosons 
[2j and Fang- Fronsdal for fermions |3]. 2 There are some approaches for the construction of 
free field Lagrangian [I] [5|. Interaction problems appear when one tries to couple massless 
higher spin fields to an electromagnetic field |6., to gravity [7j jH] or to construct self- 
interactions PO] |11| . For a review of massless higher spin gauge field theory, see |12j . 

At present, there exist various approaches to the theory, which solve the interaction 
problems in some cases. For example, an approach, called the unfolded formalism, was 
developed by Vasiliev et al |T3] [Ej. They succeeded in the construction of interacting 
higher spin gauge theory with a nonzero cosmological constant J3]. An approach, called 
BRST approach, was initiated by development of string field theory on the basis of BRST 
techniques US] [L7j [IE] H2| EH- 

In this paper, we study matrix models as a new approach to formulate massless higher 
spin gauge field theory. Recently, it has been shown that the Einstein equation can be 
obtained from the equation of motion of a matrix model by introducing a new interpretation 
of the matrix model, in which matrices represent differential operators on a curved spacetime 
PP. Furthermore, it was pointed out that there is a possibility that matrix models include 
the degrees of freedom of massless higher spin gauge fields. An advantage of this formalism 
is that the matrix model possesses a gauge invariance manifestly, embedded in the U(N) 
symmetry. Therefore it is interesting to analyze interacting massless higher spin gauge field 
theory using the matrix model. 

A first step towards constructing massless higher spin gauge field theory is the formulation 
of the free theory. Therefore, in this paper we show that the free equation of motion of bosonic 
massless higher spin gauge fields can be derived from that of the matrix model. 

There is another motivation for our study. Massless higher spin fields are expected to ap- 
pear in the tensionless limit of string theory, since mass squared of them are all proportional 
to the string tension. Thus in this limit, one should observe an enhancement of gauge sym- 

2 In four dimensional spacetime, all higher spin fields can be described either by totally symmetric tensor 
or totally symmetric tensor-spinor fields. This is not the case for a dimension that is larger than four. In 
this paper, we will restrict ourselves to the consideration of symmetric tensor fields. 



metry of string theory by that of the massless higher spin field theory. On the other hand, 
matrix models are expected to be a nonperturbative formulation of string theory. Therefore 
our study may be useful for better understanding of gauge symmetry of string theory and 
can lead to further understanding of nonperturbative aspects of string theory. 

The organization of this paper is as follows. In section 2, we briefly review the results 
of p. In section 3, we show that the free equation of motion of bosonic massless higher 
spin gauge fields can be derived from that of the matrix model. Section 4 is devoted to 
conclusions and future works. 

2 Matrix Model 

In this section, we briefly review the results of pQ. Introducing a new interpretation of a 
matrix model, we will see the following three facts. 

• Vacuum Einstein equation can be derived from the equation of motion of the matrix 
model. 

• There is a possibility that the matrix model can describe bosonic massless higher spin 
fields. 

• Gauge symmetries related to higher spin fields are embedded in the U(N) symmetry 
of the matrix model. 



In this paper, we consider the large iV reduced model of D-dimensional pure Yang-Mills 
theory with U(N) gauge symmetry as the matrix model : 

S = -^tr([A a ,A b ][A a ,A b ]) 1 (1) 

where A a are N x N hermitian bosonic matrices. Latin indices denote Euclidean spacetime 
directions. This action has the SO(D) Lorentz symmetry and N x N unitary matrix sym- 
metry. We can also consider supersymmetric version of this model, but in this paper we 
consider only the bosonic action (pQ). 3 

The basic idea of pQ is that matrices represent differential operators on a curved space- 
time. There are several problems with this identification. For example, matrices act as 
3 In D = 10, supersymmetric version of the action JIJ is nothing but the action of IIB matrix model 21 . 



Endomorphisms on a vector space, which means matrices map a vector space to itself. On 
the contrary, covariant derivatives map a tensor field of rank-n to a tensor field of rank- 
(n+1). Therefore we should prepare a vector space V which contain at least tensor fields of 
any rank and prepare an object V( a ) which is equivalent to a covariant derivative V a such 
that each component of V( a ) is expressed as an Endomorphism on V, in order to interpret 
covariant derivatives as matrices. In pQ the authors showed that V can be given by the space 
of smooth functions on the principal Spin(D) bundle over a manifold M. A smooth function 
/ on it is defined as the mapping 

f:Ux Spin(D) -> C, (2) 

where U denotes a patch on M. Thus, / depends on the local coordinate (x, g), where x G U 
and g G Spin(D), of the principal Spin(D) bundle over M. V( a ) can be given by 

V (a) = fl^V^Vb, (3) 

where R^ b (g) is the vector representation of Spin^D). 4 The covariant derivative V a is 
defined as 

V a = e^(x)(d, + u, bc (x)0 bc ), (4) 

where 0fe c (= —0 c b) is the generator of the local Lorentz group Spin(D), e a ^(x) is the vierbein 
and u}^ t ab (x) is the spin-connection. Here, Latin indices denote the local Lorentz indices. 
Notice that V a maps a rank-n tensor to a rank-(n + 1) tensor and @ a b acts on the local 
Lorentz indices of these tensors. Therefore we have 

[0a6,V c ] = -((5 ac V 6 -4cV a ), (5) 

in this setting, which will be used later. 

Let us see how the Einstein equation can be derived from the equation of motion of the 
matrix model by applying this interpretation to the matrix model. From the action (JTj), we 
obtain the following equation of motion : 

[A a ,[A a ,A b ]] = 0. (6) 



4 Though R a b (g) and R^ b (g) are the same quantity, we distinguish them because indices a is transformed 
by the action of G, while (a) is not. 



We substitute A a = iV a into the equation of motion (JSjl. where i is introduced to make A a 
hermitian. 5 The commutator of A a becomes 

[A a ,A b ] = [ieSV^iefV,,] 

= -eSefWv V„] - e [a ^V,e b] u )V u 

= —Rab C &cd — C^abVu, (7) 

where R a b,cd is the Riemann tensor and C v ' ab is the torsion tensor. Here, we put the torsion- 
less constraint C u ab = to solve the spin-connection in terms of the vierbein. Under the 
torsionless constraint, the equation of motion becomes 

[A a ,[A a ,A b ]} = \iV a ,-R ab cd O cd ] 

= -i[V\R ab cd ]O cd -iR ab cd [V\O cd ] 

= -i(V a R ab cd )O cd + iR b c V c = 0, (8) 

and we obtain the following two equations : 

R bc = 0, V a R ab , cd = 0. (9) 

The first equation is the vacuum Einstein equation. The second equation can be derived 
from the first one and the Bianchi identity : V a R bc de + V b R ca de + V ' c R ab de = 0. Thus, we 
have obtained the Einstein equation from the equation of motion of the matrix model. 

Since each component of matrices A a acts functions on the principal SO(D) bundle over 
M as an Endomorphism, in general, A a may be expanded as 

•2 

A a = iV a + a a (x) + l -{b a \x),V b } + ~{B a bc (x), bc } + l -{e a bc (x),V b V c } + ■■■, (10) 

where i and anticommutator {} are introduced to make A a hermitian. The coefficient e a bc (x) 
can be taken to be symmetric under exchange of the indices b <-> c because antisymmetric 
part can be absorbed in the term that is the first order in ab . Higher order terms expanded 
in terms of the operators V a and ab also can be taken to be symmetric under permutations 
of the operators. We consider the expansion as a sum of homogeneous polynomials in V a 
and a b, whose coefficients are identified with massless higher spin gauge fields. Coefficients 
of the first order homogeneous polynomial express spin-2 gauge fields, and those of the 



5 Hermiticity can be confirmed by denning the inner product (u,v) = j ed d xdgu*(x,g)v(x,g), where 
e = det^ea 11 ) , dg is a Haar measure of group G, u{x, g) and v(x, g) are arbitrary functions of x and g G G 



second order one express spin-3 gauge fields and so on. Thus, the number of independent 
components of higher spin gauge fields grows rapidly with degree in V a and a b- 

Here, we mention how gauge symmetries are embedded in the U(N) symmetry of the 
matrix model. Originally, the U(N) symmetry of the matrix model is written as 

5A a = i[A,A a ], (11) 

where A is a iV x N hermitian matrix. In the new interpretation of the matrix model, A 
becomes a scalar operator expanded in terms of V a and & a b- 

Let us see how gauge transformations are generated by A in the case of spin-3. In order 
to deal with this case, we need to keep track of only the first and the fifth term in (jlOj). 

A a = *V a + (*) 2 {e a , aia2 (x), V ai V a J + • • • . (12) 

We take A as 

A = (*) 2 {A aia2 (x),V ai V a2 }. (13) 

Then (fTT|) becomes 

8A a = i[A,A a ] 

= (V a A aia2 (x))V ai V a2 + --- . (14) 

Thus e a ^ aia ' 2 (x) transforms as 

6e a a ^(x) = V a X a ^(x) + --- . (15) 

Therefore, e ay aia2 (x) transforms as a rank-3 symmetric tensor field. Other gauge transfor- 
mations related to higher spin gauge fields are realized in terms of other A. Thus, this 
formulation possesses the gauge invariance related to higher spin gauge fields manifestly. 
This is an advantage of this formulation. 

3 Free Higher Spin Field Equation in Flat Spacetime 

The free equation of motion for the rank-s totally symmetric tensor field A11 - Ms in D- 
dimensional flat spacetime, the so-called "Fronsdal equation" , is given by 

w M1 ... Ms = □0 m ... Ms - sd( M i( 9 -<£U-/V) + s ( s ~ i^aA^s-^) = °> ( 16 ) 
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where we use notations (d-cp)^...^^ = d p (p p fll ... fls _ 1 and 0' M1 ... Ms _ 2 = 4> p P ^-^ s -2- Greek indices 
run from 1 to D and denote flat spacetime directions. This equation of motion possesses the 
gauge symmetry, the so-called Fronsdal symmetry, 

<%i- Ms = d^X^...^), (17) 

where the bracket () denotes symmetrization of the flat spacetime indices and the gauge 
parameter A M1 ... Ms _ 1 is symmetric under permutations of the indices. The conventional for- 
mulation for free totally symmetric tensor gauge fields was originally derived by Fronsdal 
[2|. The key feature of this formulation is the need for a pair of constraints, one on the 
parameter A Atl ... M3 _ 1 , whose trace X p PPl ...p a _ 3 = A' ..„ _ 3 is required to vanish, and one on the 
gauge field itself, whose double trace 0% CT <j Ml ... Ms _ 5 = 4>'L...u a _ B is also required to vanish. 

In this section, we show that the free equation of motion of higher spin gauge fields (jTBjl 
can be derived from that of the matrix model ©• In our formulation, as we will see in the 
next subsection, the constraints on M1 ... Ms are achieved by putting traceless constraints on 
the fields in A a . Through the analysis in this section, we can understand how higher spin 
gauge fields are included in the matrix model. 

We show the case of spin-3 in subsection 3.1 and spin-s in subsection 3.2. 

3.1 Free Spin-3 Field Equation in Flat Spacetime 

Let us first consider the spin-3 case as an example. The equation of motion of spin-3 field 
<pp Vp is given by 



W pup = D(f)p Vp -d p (d- 4>) vp -d v {d- 4>) PIX - d p {d 



/flU 



+d tl d v (/>' p + d v d P 4>' p + d p dp4>' v = 0. (18) 

In this and the next subsection, in order to derive free field equations we keep only terms 
linear with respect to the component fields and we use the ordinary differential operator <9 M 
instead of the covariant derivative V Q . Notice that the commutation relation of the operators 
is given by 

[0^,<9 P ] = ~(<W^-^A)- (19) 



In order to deal with the spin-3 case, we keep track of the second order homogeneous poly- 
nomial of the operators d p and pu in A p : 

A p = id p + i 2 e^ p {x)d v d p + i 2 ~u^°(x) {d v pa + pa d v ) 

+l ^n^( x ) (0 pa XK + XK pa ) . (20) 

Here, we do not write the anticommutator appeared in (jlOjl explicitly because terms like 
(d v eif p )d p are not needed to derive spin-3 field equation. 

Based on the analogy of the frame formulation of gravity, we can regard that the fields 
as e p>vp and u>n, u ,pa- are generalizations of the vierbein and the spin-connection, respectively 
We assume that the fields e^ up and u p ^ pa satisfy the traceless conditions 

H/ = 0, u^ p p v = 0. (21) 

The totally symmetric tensor field 4>p Vp is defined in terms of e^ V p as 

and 0Ofj, )U>P ff and VL^ pa ^ K are auxiliary fields. As we will see shortly, the relation (|2*2j) can be 
understood from the gauge transformations which are embedded in the U(N) symmetry of 
the matrix model. 

Let us summarize the gauge transformations. There are three kinds of gauge transfor- 
mations in the case of spin-3 : 

(i) Generalized coordinate transformation, generated by A = \ pv d p d v , 

OGp,vp = O^Avp, OOJft^pcr = 0, oVln^a^K = 0, (23) 

where the parameter A Ml/ satisfies A M!/ = \ up . 

(ii) Generalized local Lorentz transformation, generated by A = \\ u,pcr (d u pa + pa d v \ 

where the parameter \^ up satisfies \^ up = ~^p, P u- 
(iii) Auxiliary gauge transformation, generated by A = ^X pa ' XK (0 pa 0\K + 0\n0 P a), 

where the parameter A^ ;PCT satisfies X^pa = ~Kp,pa = -A^ ;CTp . 



Based on the analogy of the frame formulation of gravity, the gauge transformation 
(i), which corresponds to the Fronsdal symmetry, is an extension of generalized coordinate 
transformation. The gauge transformation (ii), which removes the part that is not totally 
symmetric in the indices of e PyUp , is an extension of local Lorentz transformation. The gauge 
symmetry (iii), which does not act on dynamical fields, appears only for spins larger than 2. 

Under the gauge transformation (i), the rank-3 totally symmetric tensor field (p^p defined 
in (}22j) transforms as follows : 

<5<?W = d p \ up + d v \ P p + d p \ pu . (26) 

This is consistent with the identification ()22|). 

Since the gauge transformation (i) corresponds to the Fronsdal one, in order to derive 
higher spin field equations we should fix the gauge symmetries (ii) and (iii). However, gauge 
symmetries cannot remove all extra degrees of freedom. Therefore, we must impose some 
constraints that can be used to determine the auxiliary fields 00n, v ,pcr and Qp,pa,\K in terms of 
the dynamical field e PtVp . In the rest of this subsection, we will perform these procedures in 
order. 

Constraints 

Based on the analogy of the frame formulation of gravity, we impose constraints on 
the field strengths in order to solve the auxiliary the fields LUp,u, P a and Vt^ p(J)XK in terms of 
e^ up . The field strengths are coefficients of the operators in the commutator of A p . The 
commutator is calculated as follows : 

[A„ A v \ = -id { pe v] e°d p d a - l -d { ,u u] ^ x (d p aX + aX d p ) 

- l -n [u ^ XK [d { p, {0 pu XK + O Xfi pa )] 



-f),.n., f^ XK 
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d [p n i/] ^ XK (0 p(7 XK + O XK 



pa j 



-% P ' CTX [d IM] ,dp0 aX + aX d p ] 

= -lidipe^ + u^^dpd* 

- l -(d [fl u u] ^ x + n M p ' aX + n { pr x 'u] p ){dp0,x + 0,xd p ) 
- l -d [p n u] ^ XK (0 pa XK + AK 0p,), (27) 

where the bracket [] denotes antisymmetrization of indices. We impose constraints that the 



coefficients of d 2 and 2 are equal to : 






From (|2Sjl. 



a; 



^,i/,pff is solved in terms of the first order derivatives in e p>up 



OJ 



/J,,U,p<7 



Op{ C Pt i/<r T" &tr,vfj,) 

~~Z(j(j\pp,VlL + e fl,Up) 

~T~zOp\&o,pv &p,V(j)' 

This is similar to the torsionless constraint of gravity. The constraint on Q.^pa^ 
that fiupo-AK can be written as a "pure gauge" configuration, 6 



n 



p,,pcr,\K 



C'pXpO'jXKl 



where the parameter x P <rM satisfies XpoM = ~X<t P m 
Imposing these constraints, we obtain 



~Xpcr,K\- 



1/4 : [-^/j j At 



~ [d»(d [IM uj u] ^ x + d { , x ,r x + W) 

x {d p aX + aX d p ) . 
Therefore, we can obtain the following equations of motion : 

d [p u u] ^ + d { ,Xv] M>M + d [pX ^ p) u] = 0, 

d»(d [p u y] ^ x + d [p x»] p ' aX + W) = o, 

where we define u as 



OJ 



u> 



p,(p,a)u — J^p,,p,av + ^ p.,ff,pv ■ 



P ,v,pcr — ^ P ,(p,<T)v r . VJ p,,p,(7V ■ _ - 



(28) 
(29) 



(30) 
implies 

(31) 



(32) 



(33) 
(34) 



(35) 



The equation (J51j) follows from (|3*3*|) . Therefore, dynamical field equation of motion is ()33|) . 
In order to make the equation (|53|l to be symmetric under permutations of indices u, p, o we 
should impose the following constraint on x P u, P a, 



X P u. 



pry 



o^IH.PO"- 



(36) 



'Precisely speaking, this is not pure gauge because u Pt y tPC7 is also transformed under (iii). 



Imposing this constraint, we find that the equation ()33|) is symmetric under permutations 
of indices and is second order derivatives in e^ vp . 

Here, we summarize the constraints we imposed in this subsection as follows : 

• Traceless constraints : 

e„/ = 0, u^ p % = 0. (37) 

• Field strength constraints : 

<Vh/ ,<t + ^V , *%]= ' (38) 

d\i$lv] tl »M = °' ( 39 ) 

where the constraint (JHSJ) can be solved as : fi M>pc r,AK = dpXp^M- 

• Constraint on Xp.v,po '■ 

Xpv,pe = -gW[ H , p(7 . (40) 

Imposing these constraints, the auxiliary fields w^^po an d ^Vp^Ak are expressed in terms of 
the dynamical field e^ vp and only a symmetric part remains in the equation of motion ()33j) . 
The constraints imposed on the spin-3 fields are the traceless constraints ()37|) and the 
field strength constraints (|38|). (|39|) and (|4(J|) . These constraints have been imposed in order 
to express the equation (J313|) in terms of the dynamical field e MliP2 ... Ps . However, viewed from 
the matrix model, field strengths should be introduced as independent degrees of freedom. 
There is a possibility that "the higher spin field strengths" propagate as asymmetric tensor 
fields. 

Gauge fixing 

So far, we analyzed the elimination of the extra degrees of freedom by imposing the 
constraints (|3*7jl . (|3*Sjl . (|3*9j) and (jJOJ). Combining these constraints and the equation of 
motion (J3*3*j) we find that the equation (|33J) is expressed in terms of second order derivatives 
in e p,up and is symmetric under permutations of indices. However, these constraints cannot 
eliminate all extra degrees of freedom. The Fronsdal equation (|18|) is expressed in terms of 
the rank-3 totally symmetric tensor field 4> pupi but the equation (|3*3*|l is expressed in terms of 
"the spin-3 vierbein" e^up, which have the part that is not totally symmetric. Last remaining 

10 



extra degrees of freedom is the part that is not totally symmetric in the indices of e^ up . Thus, 
we should eliminate this extra degrees of freedom and express dynamical variable in terms 
of (f)fj,u P , in order to show that the equation (jHSJ) coincides with the Fronsdal equation (JTBJ). 
This can be done by fixing the gauge symmetries. Recall that there are three kinds of the 
gauge transformations (i), (ii) and (iii) and the gauge transformation (i) corresponds to the 
Fronsdal gauge transformation. Therefore it seems that we should fix the gauge symmetries 
(ii) and (iii). As we see, the gauge symmetries (ii) and (iii) can eliminate the part that is 
not totally symmetric in the indices of e P)Up and we can express dynamical variable in terms 
of the rank-3 totally symmetric tensor field : (p^p. 

First, we fix the gauge symmetry (iii). Gauge fixing can be done by transforming cD^pcr ~~ > 
^n,v,pa = &n,v,po + ^^,pa, choosing the parameter A F|/W as 

**-/jU/,p<T = ^fi,u,pa ~r ~ -Dpi/, pa ~7\^pp,av t LJpo,pv *5pv,po ^av,pp) 

-PVpp(B'^ u + B'^) - (3 VfMT (B ptl/ + BlJ 
-PVpu(B'^ + B'^) - ^ v {B'p,p + B^) 

+2P V , U (B'^ + B' fftp ) + 2/V(ai,* + S i,,)> ( 41 ) 

where we define B pvpa , B' and (3 as follows : 

*5pv,po = ^°Jp,u,p<j ^u,p,pa 

(Jp&u,p<T (Jv&p,pcri V^^l 

B'„ = B„%, (43) 

" = W=Ty (44 > 

^pv,po satisfies the following properties : 

l*-pv,p<j l*-vp,p<j l*-pii,tjpi ^p{v,pa) U, x*^) 

Kp% = °> V/ = °- ( 46 ) 

(3 is determined by the traceless condition ([46)1 . 

Next, we fix the gauge symmetry (ii) by transforming e M)1/p — > e^ up = e p ^ p + A PiUp with 

Ap,u P = -e p ,u P + <\>pv P + a(Vp»<P'p + Vpp<P'u - 2 Vvp<l>'p)- (47) 

Ap,u P satisfies the following properties : 

Ap,up — ^p,pv, ^(p,vp) — 0, (48) 

K/ = 0, (49) 

11 



where 0^ = 4> pp p , a = 2 (d-x) • a * s determined by the traceless condition (jUJ). Carrying out 
these transformation, we can remove the part that is not totally symmetric in the indices of 
e^ up and we have e M)l/p = 4>^ up . Substituting w^ v ^ pa and e^p into the equation (JHBJ), we can 
show that the equation (f3"3j) coincides with the Fronsdal equation (JTHJ). 

It is worth noting the relation between our formulation and the unfolded formalism of 
higher spin gauge field theory due to Vasiliev [13] • Vasiliev constructed free field Lagrangian 
using similar method to the one we have employed in this paper [Hj. Only difference between 
the unfolded formalism and our formulation is the appearance of O^^a* i n free theory. In 
the unfolded formalism, ^l^ pa% \ R appears in interacting theories and contribute to higher 
derivative interactions. It is interesting to investigate the relation between our formulation 
and the unfolded formalism by analyzing higher spin interactions. 

3.2 Free Spin-s Equation in Flat Spacetime 

In this subsection, using the same method as the one we have employed in the previous 
subsection, we derive the free equation of motion of the rank-s totally symmetric tensor field 
in D- dimensional flat spacetime (|T6|) from that of the matrix model. 

In order to deal with the spin-s case, we keep track of the (s — l)-th order homogeneous 
polynomial of the operators d p and Mi , : 

* _ -I "7*1 \ u Pl u Pa-2 VJ P\d\\ 

ay- 1 

+ l2 O t*i"Hs-3,pi<n,p2<T2fn ...a i 

( — 1V — 9 1 ^' \ U Pl U Ps-'i VJ PlOl VJ P2<72} 

ay- 1 

4- 12 O.,. P-l-Ps-i,Pl<Tl,P2<r2,P3(T3rfl ...f) 1 

( s _ 1 Ve _ OVs — 3) (y Ms-4 V ^Pl0-i^'p20-2 Vy P30-3j 

+ 77^^(-3)/ 1<ri '"''' Ps - 1<Ts - 1 {0^ 1 ---0p s -^_ 1 }. (50) 

^Ati-zxa-i^) is the "spin-s vierbein" and w M , M1 ... /Xs _ 2 ,po-(a ; ) is "the spin-s connection". From 
the discussion in the previous subsection, it seems that Cl^(i = 1, • • • , s— 3) are not necessary 
to derive the equation of motion. We set these auxiliary fields to zero : Qi = ■ ■ ■ = f2 s _3 = 0. 
Therefore, A^ becomes 
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A, = id, + (i)-V' 1 ~'*'- 1 M ...0 M ._ 1 



ay- 1 



(i) 8 ' 1 

(c _ 1 V s _ 9) A *' Pl a M s -3^Pl0-l^P20- 2 J- \° l ) 

Here we assume that the spin-s fields satisfy the traceless condition, 

The rank-s symmetric tensor field </> pl ... Ps is defined as M1 ... Ms = e( AtliAta .„ /Xa ). 0^...^ satisfies 
the double traceless condition p p CT o - M5 ... Ais = as a consequence of the traceless condition 

452). 

Let us summarize the gauge transformations. In the case of spin-s, there are s kinds of 
gauge transformation. Since auxiliary fields fl(i)(i = 1, • • • , s — 3) are set to be zero, the 
following three kinds of gauge transformations remain : 

(i) Generalized coordinate transformation, generated by A = A M1 *" Ms_1 9 M1 • • • Ms _ 1 , 

O e /j,/Lti— Ms-l = 0'/xA w ...^ s _ 1 , OWfi,/i I -/i s _ 2 ,/jj_ii' = U, oilij itl j, 1 ... f j, a _ 3t p (Tt x K = U, (ooj 

where A^...^^ is symmetric under permutations of the indices /ii, • • • , fJ, s -i- 

(ii) Generalized local Lorentz transformation, generated by A = A Ml '" Ms - 2,po "{<9 Ml • • • <9 Ms _ 2 pCT }, 

" e M.Mi— Ms-i ~~ / V^i--/x s _2,Ais-i)M = ''ViMl—Ms-1' 

S£lwv"f.-3,peM = °> ( 54 ) 

where A /U1 ... /Us _ 2i p er is symmetric under permutations of the indices /Ji, • • • , /i s -2 and 
saiisnes A pl ... Ps2 ,pcr A P i--- Ps _2,<rp- 

(iii) Auxiliary gauge transformation, generated by A = A jUl ... jUs _ 3i/9O . : A K {0 Atl • • • 9 Ms _ 3 pct Ak }, 

" e p,Pi--p s _i = u, 

"^M.Ml— Ms-2,P°" ~~ / ^(/ll---/l 3 _3,p s _2)p,po- + ^(yUl— |Us-3,PC,Ms-2)M> 

0"/i,|*i— /j s _ 3 ,pcr,Are = 6> / Vi-Ms-3,P°",Ak; I' 3 ' 3 ) 

where Xfj, v ..fj, a _ 3 ,per,XK is symmetric under permutations of the indices /xi, • • • , /i s -3 and 
satisfies A Pi ... Ps3i/9CTj ak = — A Pl ... Ps3jCr p,AK = ~~ A^ 1 ... /is _ 3i p (T) «A- 
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Therefore 



) YfJ,i---fJ. B 



e (/*Li/*2"7»») 



are transformed under (i) as 



'A*!— Ms 



^(Ml ^M2—Ms) • 



(56) 



This is consistent with the identification 0^...^ = e^^...^). In order to derive the equation 
of motion (|16|) from that of the matrix model, what we must do is to impose constraints on 
the auxiliary fields and fix the gauge symmetries (ii) and (iii), as in the spin-3 case. 

Constraints 



Constraints are imposed on the coefficients of d s and d s in the commutator of Ap. 
We obtain 






(57) 
(58) 



The constraint (J55j) implies that ^l lh fi 1 ... lls _^ pa> \ K can be written as a "pure gauge" configu- 
ration : 



n 



/l,tll---/ls-3,pO-,XK 



°>X/ii— 



fJ, s -3,p(T,\K,- 



(59) 



We impose constraints on Xfiv-fi a -z,paM as 



Xa*i--m s -3,PO',Ak 



-a; 



[p,o-]^ r --/i s _3,AK- 



(60) 



Imposing these constraints, the commutators of A^ become 



[Afj, , Ay\ 



8 ~ 1 



a M ^i^-^- 2 '^ - a/i"^-3^-^ ^ - n 



v 



[*« 



/il---/i s _3,po-,/i s _2 



xRr-^0,,}, (61) 



and 



[^MA^A,]] 



3[ /i u; I/ ]' ir "' i - a ' M - 1 ' i 



Q /il---// s -3,A i s-2 Ms-lM _ Or A i l'"Ms-3,Ms-lM ) Ms-2 



xRi---^ s -i} 



^(^w^" 1 -"- 8 '^ 



-fil 



Ml---A i s-3,A t s-2 , PCT 



a^»^^^ x {9 . . .3 }.( 6 2) 
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Therefore we obtain the equations of motion 

dyjj^W"*-* 11 - 1 - Q { ^--n s - 3 ,ns-2^ns-i)n _ fi ^(w-M.-3,/i.-i/i,/i.-3) j/] = 0) (63) 

d" (d^oj^ 1 -^- 2 '^ - n { / 1 '"»- 3 >' l -\ ]t '" - n^ 1 "^- 3 '^^- 9 ,,]) = 0, (64) 

where Cj^^...^^ is defined as 

The second equation ()64j) can be derived by using the first one (|63|). Therefore, dynamical 
field equation is the first one. Owing to the constraint (|60|). the equation of motion ([63)1 are 
symmetric under permutations of indices. 

Here, we summarize the constraints we imposed in this subsection as follows : 

• Traceless constraints : 

• Field strength constraints : 

V"]." 1 '""' -1 + ^r"^- 2 ^' 1 ^] = 0, (67) 

^X^I/] ) Atl-/i._3,P<T,AK = 0. (68) 

The constraint ((58)) can be solved as 

• Constraints on Xm-m s _ 3 ,p<t,Ak : 

Imposing these constraints, the auxiliary fields u) lii n 1 ... fia _ 3tPff and fi jtt ,^ 1 ... A4B _ 3 ,^ T) AK are expressed 
in terms of the dynamical field e MjM1 ... Ms _ 1 and only a symmetric part remains in the equation 
of motion (jHSJ). 

Gauge fixing 
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First, we fix the gauge symmetry (iii) by transforming u^^...^^ -> w^^...^^ = 
Vu,n,ttvi*.-i +^M.Mi-Ma-i> choosing the parameter Kw-it—i 

+ Pl^ A t|-D A j liAt2 ... Msl |/i i 

+ P21 ? ?^/ii-t> A t )At2 ... / _ ts _ 1 + r ]^^l^u,fl2---fJ.s-lJf J 'i 

— {Pi + P2j|^/xi-0 ;U2)W 3... A i s _ 1 + ? ?wn-"M2,^3---Ms-iJ>> 

+(A + 2/3 2 ){r ?w/i2 S; 3iW4 ... /Xs _ 1 } w , (71) 

where {} w denotes the symmetrization of indices /i«(i = 1, • • • , s — 1). We define S AM/)Atl ... /Xs _ 1 , 
S mi,M2-m s -i' & and ^ 2 as follows : 

— o fl e u ^ 1 ... fls _ 1 — o u e fl ^ 1 ... fls _ 1 , {(2) 

A - ^|, (74) 

^ 2 " 2 (D-2)(3- S -D)- (75) 

The parameter A VMjM1 ...^ j) _ 1 is symmetric under permutations of the indices fii,- • • , fJ, s -i and 
satisfies the following properties : 

A p = o a p — n C77 , i 

/?i and /?2 are determined by the traceless conditions (|77|). 

Next, we fix the gauge symmetry (ii) by transforming e MliM2 ... Ms -> £ M1)M2 ... Ma = e Mll/J2 ... ft + 
A WlW ... M . with 

where a = 2(1+ J-4) • ^-ihw-ij* * s symmetric under permutations of indices /i 2 ,- • ■ , /i s and 
satisfies the traceless condition 
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Carrying out these transformation, we can remove the part that is not totally symmetric 
in the indices of e^^...^ and we have e WiW ... fe = (p^-^- Substituting w^^...^^ and 
£ m,tj.2-na m ^° the equation (pHJ) . we can show that (jHSJ) coincides with the Fronsdal equation 

(HE). 

4 Conclusions and Future Works 

In this paper, we have shown that the free equation of motion of bosonic massless higher spin 
gauge fields in D-dimensional flat spacetime can be derived from that of the matrix model 
based on the new interpretation of the matrix model. In order to derive higher spin field 
equations, we have done the two things : 1) imposing constraints 2) performing gauge fixing 
procedures. The results of this paper suggest that bosonic massless higher spin fields can 
be embedded in the matrix model. This is a first step towards construction of interacting 
massless higher spin gauge field theory by using the matrix model. Therefore, the matrix 
model can be used as a new approach to formulate massless higher spin gauge field theory. 
There are several things which should be studied further. One is the derivation of free 
fermionic massless higher spin gauge field equations. Recently, it was shown that supergrav- 
ity can be embedded in the supermatrix model [22] • There is a possibility that fermionic 
higher spin fields are embedded in the supermatrix model. Another one is to construct the 
interacting massless higher spin gauge field theory. As mentioned in Introduction, it is dif- 
ficult to construct interacting theory. Difficulties associated with the requirement of gauge 
invariance can be overcome by using the matrix model because it has gauge invariance man- 
ifestly. 
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